Long range scattering for the Maxwell-Schrodinger system 
with arbitrarily large asymptotic data 
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Abstract 



We review the proof of existence and uniqueness of solutions of the Maxwell- 
Schrodinger system in a neighborhood of infinity in time, with prescribed 
asymptotic behaviour defined in terms of asymptotic data, without any re- 
^ | striction on the size of those data. That result is the basic step in the con- 

struction of modified wave operators for the Maxwell-Schrodinger system. 
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1 Introduction 



This paper is devoted to the theory of scattering and more precisely to the 
construction of modified wave operators for the Maxwell-Schrodinger (MS) system 
in space dimension 3, namely 

' id t u = -(1/2)A A « + A e u 

< DA e -d t (d t A e + V-A) = \u\ 2 (1.1) 
OA + V {d t A e + V • A) = Im uS7 A u 

where u and (A, A e ) are respectively a complex valued function and an IR 3+1 valued 
function defined in space time IR 3+1 , = V — %A and A^ = are the covariant 
gradient and covariant Laplacian respectively, and □ = d 2 — A is the d'Alembertian. 
An important property of that system is its gauge invariance, namely the invariance 
under the transformation 

(it, A, A e ) -> (it exp(-i0), A - V6>, A e + d t d) , 

where 9 is an arbitrary real function defined in IR 3+1 . As a consequence of that 
invariance, the system (11.11) is underdetermined as an evolution system and has to 
be supplemented by an additional equation, called a gauge condition. Here, we shall 
use exclusively the Coulomb gauge condition, V • A = 0. Under that condition, the 
equation for A e can be solved by 

A e = -A _1 |?i| 2 = (47r|a;|)" 1 * \u\ 2 = g(u) (1.2) 

where * denotes the convolution in IR 3 . Substituting ( 11.2j) and the gauge condition 
into ( 11. II) yields the formally equivalent system 

' id t u = -(1/2) A A u + g(u)u (1.3) 

< 

□A = P Im «Va« (1.4) 

where P = 1 — VA _1 V is the projector on divergence free vector fields. 

The MS system is known to be locally well posed both in the Coulomb gauge 
and in the Lorentz gauge d t A e + V • A = in sufficiently regular spaces [7] [8] , to 
have weak global solutions in the energy space [6] and to be globally well posed in 
a space smaller than the energy space [9]. 
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A large amount of work has been devoted to the theory of scattering and more 
precisely to the existence of wave operators for nonlinear equations and systems 
centering on the Schrodinger equation, in particular for the Maxwell-Schrodinger 
system PQ [2] [5] [HI] [II]- As in the case of the linear Schrodinger equation, 
one must distinguish the short range case from the long range case. In the former 
case, ordinary wave operators are expected and in a number of cases proved to 
exist, describing solutions where the Schrodinger function behaves asymptotically 
like a solution of the free Schrodinger equation. In the latter case, ordinary wave 
operators do not exist and have to be replaced by modified wave operators including 
an additional phase in the asymptotic behaviour of the Schrodinger function. In that 
respect, the MS system in IR 3+1 belongs to the borderline (Coulomb) long range case. 
General background and additional references on that matter can be found in [3] 



The main step in the construction of the (modified) wave operators consists in 
solving the local Cauchy problem with infinite initial time and we shall concentrate 
on that step. In the long range case where that problem is singular, it amounts to 
constructing solutions with prescribed (singular) asymptotic behaviour in time. The 
results include the uniqueness and the existence of solutions with such prescribed 
asymptotic behaviour. The construction is performed by changing variables from 
(u, A) to new variables, replacing the original system by an auxiliary system for the 
new variables, solving the corresponding problems for the auxiliary system and fi- 
nally returning therefrom to the original one. The auxiliary system will be described 
in Section 2 below, and the uniqueness and existence results will be presented in 
Sections 3 and 4 respectively. The exposition is based on [3] [5]. 

2 The auxiliary system 

In this section we derive the auxiliary system that will replace the original system 
(1.3) (1.4). We first replace the equation (1.4) by the associated integral equation 
with prescribed asymptotic data (A + ,A + ), namely 



0- 




Jt 



where u 



(-A) 



/ 2 , Aq is the solution of the free wave equation DAq = given by 




(2.2) 
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and we assume V -A + = V- A + = in order to ensure the gauge condition V-A = 0. 
We then perform a change of unknown functions which is well adapted to the study 
of the system (1.3) (12 .ip for large time. The unitary group U(t) which solves the 
free Schrodinger equation can be written as 

U(t) = exp(z(*/2)A) = M(t) D(t) F M(t) 

where M(t) is the operator of multiplication by the function 

M{t) = exp (ix 2 /2t) , 

F is the Fourier transform and D(t) is the dilation operator defined by 

D(t) = (it)- 3 ' 2 D (t) , (D (t)f) (x) = f(x/t) . 

We first change u to its pseudo conformal inverse u c defined by 



u(t) = M(t) D(t) u c (l/t) (2.3) 

or equivalently 



u(t) = Fu c (l/t) , (2.4) 
where for any function / of space time 

f(t,-) = U(-t) /(*,-) . (2.5) 

Correspondingly we change A to B defined by 

B{t) = -t- l D Q {t) A(l/t) . (2.6) 

The transformation (u, A) — > (u c , B) is involutive. Furthermore it replaces the 
study of [u, A) in a neighborhood of infinity in time by the study of (u c , B) in a 
neighborhood of t = 0. 

Substituting ( 12. 3ft ( 12.61) into (1.3) and commuting the Schrodinger operator with 
MD, we obtain 

{(id t + (l/2)A A -g(u))u}(t) 

= r 2 M(t) D(t) {{id v + (l/2)A B(t0 - B(t') - t'- l g{u c {t')))u c {t')} ti=i/t 
where for any IB? vector valued function / of space time 

f(t,x) = r 1 x-f(t,x) . (2.7) 
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Furthermore 

Im («Vi«) (t) = r 3 D (t) {x\u c {t')\ 2 - t' (Im u c V B u c ) (*')} t , =1 

by a direct computation, so that the system (1.3) (12. ip becomes 
id t u c = -(l/2)A B u c + Bu c + i -1 c?(u c )u c 

B 2 = B 2 (u c ,B) 

where 



A 



(2.8) 
(2.9) 



B 2 = B - B - B l , (2.10) 

B (t) = -r x D Q (t) A Q (l/t) , (2.6) 

B 1 ^B 1 (u c ) = -F 1 (Px\u c \ 2 ) , (2.11) 

B 2 (u c , B)=tF 2 (P Im u c V B u c ) , (2.12) 

/oo 
z/ -2 ^ uj- 1 sm(u(v - l))D (v) M(t/u) . (2.13) 

Here we take the point of view that B\ is an explicit function of u c defined by (2.11) 
and that (2.10) is a change of dynamical variable from B to B 2 . The equation (2.9) 
then replaces (12.11) . 

In order to take into account the long range character of the MS system, we 
parametrize u c in terms of a complex amplitude v and a real phase ip by 



Uc 



(2.14) 



The role of the phase is to cancel the long range terms in (2.8), namely the contri- 
bution of B\ to B and the term t _1 5f(M c ). Because of the limited regularity of B\, 
it is convenient to split B\ and B into a short range and a long range part. Let 
X G C°°(IR 3 , IR),0<x<l, X(0 = 1 for If I < 1, x(0 = for \£\ > 2. We define 

B L = B 1L = F* X (- t^F B x 

(2-15) 

Bs = B + Bis + B 2 , Bis = Bi — B 1L , 



We then obtain the following system for (v, ip, B 2 ) 

idfV = Hv 

dw = t- 1 g{v) + Bi L {v) 
{ B 2 = B 2 (v,K) 
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(2.16) 
(2.17) 
(2.18) 



where 

H = -(1/2)A K + B S , (2.19) 

K = B + Vip, (2.20) 

by imposing (2.17) as the equation for <p. Under (2.17), the equation (2.8) becomes 
(2.16). The system (2.16)-(2.18) is the auxiliary system which replaces the original 
system (1.3) (1.4). 

3 Uniqueness of solutions 

In this section we present the uniqueness result at infinity in time for the MS 
system in the form (1.3) (12.11) . Since the Cauchy problem for that system is singular 
at t = oo, especially as regards the function u, the uniqueness result for that system 
takes a slightly unusual form. Roughly speaking it states that two solutions (««, Aj), 
i = 1,2, coincide provided Ui and A { — A do not blow up too fast and provided 
«i — U2 tends to zero in a suitable sense as t — > oo. In particular that result does 
not make any reference to the asymptotic data for u, which should characterize its 
behaviour at infinity. 

We denote by || ■ || r the norm in U = L r (IR 3 ), by H k the standard Sobolev 
spaces and by H k their homogeneous versions, with H 1 C L 6 . We need the spaces 

V = {v : v G H 3 , xv G H 2 ) , (3.1) 
V+ = [v :< x > 3 v E L 2 , < x > 2 Vv G L 2 } (3.2) 

which are Fourier transformed of each other, and the dilation operator 

S = td t + x-V + l. (3.3) 

We denote nonnegative integers by j, k. 

We first state the uniqueness result for the auxiliary system (2.16)-(2.18). 

Proposition 3.1. Let < r < 1, I — (0, r] and a > 0. Let A be a divergence free 
solution of the free wave equation UAq = such that Bq defined by (2.6)q satisfy 

II V k (td t yB |U + || W k B |U < C r k (3.4) 

for < j + k < 1 and for all t G /. Let (t> i? (ft, B 2 i), i = 1,2, be two solutions of 
the auxiliary system (2.16)-(2.18) such that v t G L% C (I,V), Wifi G L% C (I , H 1 f] H 2 ) , 
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B 2i G L^ C (J, if 1 ), satisfying the estimates 

in t) a , (3.5) 

in t) 2a , (3.6) 

2 < C t 1+£ (3.7) 

for some e > and for all t G I, and such that (pi(t) — (p 2 {t) tends to zero in H 1 
when t — > 0. T/ien (vi, tpi, B 21 ) = (v 2 , f 2 , B 22 ). 

Remark. Although each fi separately is not in H 1 , it follows from the equation 
(2.17) and from the assumption (13.71) that tpi — tp 2 has a limit in H 1 when t — > 0, 
which gives a meaning to the assumption that (pi — (p 2 tends to zero in H l when 
t -> 0. 

We next state the uniqueness result for the MS system in the form (1.3) (12.11) . 
We recall that u(t) = U{-t)u(t). 

Proposition 3.2. Let 1 < T < oo, I = [T, oo) and a > 0. Let Aq be a divergence 
free solution of the free wave equation DAq = satisfying 

|| V k S j A (t) ||oo + || V k x ■ A (t) ||oc < C r 1 (3.8) 

for < j + k < 1 anc? /or a// t G I. Le£ (uj, Aj) ; z = 1, 2 ; be two solutions of the 
system (1.3) / TO) suc/i £/&a£ u { G L^ C (I, K), - A) G L^ c (I,H l ), satisfying the 
estimates 

II «.(*); K || < C(l + £n t) a , (3.9) 

|| V(A - A )(t) || 2 < C r 1/2 (l + £n t) 2Q , (3.10) 

||< x/t > (ut - u 2 )(t) \\ 2 < C r 1_£ (3.11) 
/or some e > and /or a// 1 G J. T/ien (wi, Al) = («2, ^2)- 

The proof of Proposition 3.1 follows from an estimate of the difference of two 
solutions of the system (2.16)-(2.18), in particular of V\ — v 2 in H 1 C\F(H 1 ). Proposi- 
tion 3.2 is proved by reducing it to Proposition 3.1. The main step of that reduction 
is the construction of the phases <pi, and that step requires a strengthening of the 
assumptions in the form of the replacement of a by 7a + 3 from (13.91) (13. 10[) to 



II Vi(ty,v\\ <c(i- 

II VB^t) \\ 2 < C(l - 

< x > (vi(t) -v 2 (t)) II 
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(13.51) (13. 6p . The assumptions (13 .4p on B and (13.81) on A are equivalent. Since SA 
and x ■ A are solutions of the free wave equations, those assumptions can be im- 
plemented under suitable well known conditions on the asymptotic data (A + ,A + ). 
Finally, Proposition 3.2 is sufficient to cover the case of the solutions of the MS 
system constructed in the next section. 

4 Existence of solutions 

In this section, we sketch the proof of existence of solutions of the MS system 
(1.3) (12. ip with prescribed asymptotic behaviour at infinity in time, and for that 
purpose we first construct solutions of the auxiliary system (2.16)-(2.18) with pre- 
scribed asymptotic behaviour at time zero. Since that system is singular at t = 0, 
we proceed by an indirect method. We pick an asymptotic form for the dynam- 
ical variables, to be chosen later, we rewrite the auxiliary system in terms of the 
difference variables and we solve the latter for those variables under the condition 
that they tend to zero at time zero. Let (v a , (p a , Bi a , B 2a ) be the asymptotic form 
of (v , if, Bi, B2) and define the difference variables 

(w,tl>,Gi,G 2 ) = (v,cp, B 1 ,B 2 ) - (v a , ip a ,B la , B 2a ) , (4.1) 

where we do not assume that = Bi{v a ) or ip a = ip(v a ), in order to allow for 
more flexibility. Substituting (14. ip into the auxiliary system (2.16)-(2.18), where in 
addition we reintroduce an equation for Bi, we obtain the following system for the 
difference variables : 

' id t w = Hw + H lVa - i?i (4.2) 

d t ip = t~ l g(w, 2v a + w) + B 1L (w, 2v a + w) - R 2 (4.3) 

< 

G 1 = B l (w,2v a + w)-R 3 (4.4) 

^ G 2 = B 2 (w, 2v a + w,K) + b 2 {v a , L) - R 4 (4.5) 

where H, K are defined by (2.19) (2.20), K a = B a + Vip a = B + B la + B 2a + V<p a , 
L = Gt + G 2 + W, B x and B 2 are defined by (2.1 1)-(2. 13) , 

b 2 (v a ,L)=B 2 (v a ,K)-B 2 (v a ,K a ) , 

so that b 2 (v a , L) is linear in L, 

H l= iL- V Ka + (i/2)AV + (1/2)L 2 + G 1S + G 2 , (4.6) 
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Gis and G 2 are defined according to (2.7) (2.15), and for any quadratic form Q(f), 
we denote by Q(fi, f 2 ) the associated polarized sesquilinear form. The remainders 
Ri, 1 < i < 4, are defined by 

Ri = id t v a - H a v a = id t v a + (l/2)Atf w a - i?aS 
#2 = <9 t y? a - t" 1 ^^) - B 1L (v a ) 



R 3 = B la -B 1 (v a ) 
[ R± = B 2a -B 2 (v a ,K a ) . 



(4.7) 
(4.8) 
(4.9) 
(4.10) 



They measure the failure of the asymptotic form to satisfy the original system (2.16)- 
(2.18) and their time decay as t — > measures its quality as an asymptotic form. 
The resolution of the new system (4.2)-(4.5) now proceeds in two steps. 



Step 1. One solves the system (4.2)-(4.5) for (w, ip, G 1: G 2 ) tending to zero as t — > 
under general boundedness properties of (v a , ip a , B\ a , B 2a ) and general decay assump- 
tions on the remainders Ri, 1 < i < 4, as t — > 0. For that purpose, one linearizes 
the system (4.2)-(4.5) partly into a system for new variables (w',G' 2 ), namely 

id t w' = Hw' + Hiv a — Ri 

(4.11) 

G' 2 = B 2 {w, 2v a + w,K) + b 2 {v a , L)-R 4 . 

(There is no point in introducing new variables for if) and G\ since they are explicitly 
given as functions of w under the condition ^(0) = 0). For given (w,G 2 ) tending 
to zero as t — > 0, one solves the system (4.11) for (w',G' 2 ) also tending to zero at 
t — > 0. This defines a map T : (w, G 2 ) — > (w', G' 2 ). One then shows by a contraction 
method that the map Y has a fixed point. 



Step 2. One constructs asymptotic functions (v a , f a , B la , B 2a ) satisfying the assump- 
tions needed for Step 1. For that purpose, one solves the auxiliary system (2.16)- 
(2.18) approximately by an iteration procedure. It turns out that the second ap- 
proximation is sufficient. 

An essential ingredient in the implementation of Step 1 is the choice of the 
function space where the map V is to have a fixed point. The definition of that 
space has to include the local regularity of the relevant functions (v, B 2 ) or (w, G 2 ) 
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and the decay of (w, G 2 ) as t — > 0. As regards the local regularity of v or w, we 
shall use again the space V defined by (3.1), namely we shall take v, w G C(J, V) for 
some interval I. On the other hand, when dealing with the Schrodinger equation, 
one time derivative is homogeneous to two space derivatives, so that a given level 
of regularity in space can be obtained by estimating lower order derivatives if time 
derivatives are used. In the present case, it turns out that regularity of v is required 
at the level of H k with k > 5/2, and this is conveniently achieved by using one space 
and one time derivative, namely with v G C 1 (/,iJ 1 ). The local regularity of B 2 or 
G 2 is less crucial and is essentially dictated by the available estimates. Summing up, 
the local regularity is conveniently encompassed in the definition of the following 
space, where / C (0, 1] : 

X (I) = {(v, B 2 ):ve C(I, V) n C\I, H 1 n FH 1 ) , 

B 2 , B 2 G C(I, H 1 n H 2 ) n C\I, H 1 )} . (4.12) 

The function space where the map Y is to have a fixed point should also include 
in its definition the time decay of (w, G 2 ) as t — > 0. For that purpose, we introduce 
a function h G C(I,IR + ), where / = (0, r] such that the function h(t) = t~ 3 / 2 h(t) 
be non decreasing and satisfy 

/' dt' t'" 1 hit') < chit) (4.13) 
Jo 

for some c > and for all t G /. Typical examples of such a function would be 
hit) = t 3//2+A (l — in t)^ for some A > 0. The function hit) will characterize the time 
decay of || w(t) \\ 2 as t — > 0. Derivatives of w will have a weaker decay, typically 
with a loss of one power of t for a time derivative, and of half a power of t for a 
space derivative. The time decay of G 2 is tailored to fit the available estimates, with 
a loss of one power of t per time or space derivative. Finally, the relevant space for 
(w, G 2 ) is defined by 

X(I) = Uw,G 2 )eX (I):\\ (w,G 2 );X(I) || = Sup hit)' 1 N(w, G 2 )(t) < oo\ 
L tei J 

(4.14) 

where 

N(w,G 2 ) = \\< x > w \\ 2 V t (||< x > d t w \\ 2 V ||< x > Aw \\ 2 ) 
V t 3/2 (|| Vd t w \\ 2 V || VAw \\ 2 ) v r 1/2 1| VG 2 || 2 
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v t 1 ' 2 (n V 2 G 2 || 2 V || VO t G 2 || 2 V || VG 2 || 2 ) 

V e' 2 (II V 2 G 2 || 2 V || Vd t G 2 || 2 ) . (4.15) 

We now turn to the implementation of Step 1 and for that purpose we need 
general boundedness properties of (v a , f a , B a ) and general decay assumptions of the 
remainders Ri, l<i<4ast^0, which we state as follows. Here / is some interval 
I = (0,r ] with < r < 1. 

(Al) Boundedness properties of v a . We assume that 

v a e (CnL°°)(I,V) (4.16) 

t l / 2 d t v a e (Cnr)(/,F 2 ) , t l ' 2 xd t v a e (C n L°°)(l, H l ) . (4.17) 



(A2) Boundedness properties of <p a , B a . We define s a = Vy? a and we recall that 
K a = B a + s a = B + B la + B 2a + s a and that B aS = B + B laS + B 2a in analogy 
with (2.20 and (2.15). We assume that the following estimates hold for all t G / : 

II K a IU < C(l-£nt) , (4.18) 

II d t K a IU V || VK a IU V * || Vd t K a IU < C r 1 , (4.19) 
V || VV • s a || 3 V || V(B la + 5 2a ) Hoc 

v f ( || vd t s a iu v || vd t v ■ s a || 3 v || vd t {B la + 5 2a ) ||oo ) < c r 1 ' 2 , (4.20) 
II VB a iu v * || vd t B a iu <cr', (4.21) 

II B aS Hoc V i || d t B aS IU < C t- 1 / 2 . (4.22) 

(A3) Time decay of the remainders. We assume that the remainders Rj satisfy 
the following estimates : 

(||< x > R^t) || 2 <) ||< x > $#1; £*(((), f],L 2 ) || < n r 1 /i(f) , (4.23) 

(II Vi2i(*) || 2 <) || Vc^L^O^L 2 ) || < n r 3/2 , (4.24) 
|| V k+1 R 2 (t) || 2 < r 2 t" 1 "^ /i(t) for k = 0, 1, 2 , (4.25) 

II vi? 3 1| 2 v t 1 ' 2 1| v 2 i? 3 lh v t (ll V9 t i? 3 || 2 v || Vi? 3 Ih) 

V t 3/2 || V 2 i? 3 || 2 V t 2 || V^i? 3 || 2 < r 3 /i(t) , (4.26) 
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II Vi? 4 || 2 V t (ll V 2 i? 4 || 2 V II Vd t R A || 2 V || Vi? 4 || 2 ) 

V t 2 ((I V 2 i? 4 lb V II Vd t R 4 || 2 ) <r 4 t 1/2 (4.27) 
for some positive constants rj, 1 < j < 4 and for all t £ I. 

The decay properties of the remainders are tailored to fit the definition of the 
space X(I). For instance when applied to (4.2), the estimate (4.23) yields a contri- 
bution 0(h(t)) to ||< x > w(t) || 2 by integration over time. 

In order to ensure that B satisfies the properties of B a appearing in (A2), we 
need some regularity assumptions on the asymptotic data (A + , A + ). We shall even- 
tually make the following assumption, which suffices for that purpose. 

(A4) Regularity of (A+,A+). We assume that V • A + = V ■ A + = and that 
(A|_, A + ) satisfies the conditions 

A G L 2 , V 2 A G L 1 , uj- 1 A G L 2 , V^i G L 1 (4.28) 

for 

r ^ = (:r • V)- ? 'V fe A+ ^ = (x ■ VyV k (x ■ A+) 

(4.29) 

[ A={x- V) j V k A + A={x- V) j V k {x ■ A+) 

for < j,k < 1. 

We can now state the main result of Step 1. 

Proposition 4.1. Let the assumptions (Al)-(AS) be satisfied. Then there exists t, 
< r < 1, and there exists a unique solution (w,ip,G2) of the system (4-2)-(4-5) 
such that (w, G2) G X((0,r]) and that ip(0) = 0. Equivalently, there exists a unique 
solution (v, (p, B 2 ) of the system (2.16)-(2.18) such that (v—v a , B 2 — B 2a ) G A((0, r]) 
and such that (ip — <p a )(t) ^0 as t — > 0. 

We now turn to Step 2, namely the construction of (v a , ip a , B a ) satisfying the 
assumptions (A1)-(A3). Taking for orientation h(t) = t 3 / 2+A with A > 0, we need in 
particular that 

|| R^t) h = 0(t 1 ^) 

in order to ensure (4.23). The obvious choice v a = U(t)v + yields y? a = 0(1 — in t) 
when substituted into (2.17), and therefore Ri = 0((1 — in i) 2 ) when substituted 
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into (4.7), which is off by slightly more than half a power of t. We need therefore 
to go to the next approximation in the resolution of the auxiliary system (2.16)- 
(2.18). Omitting for the moment the terms containing B , which require a separate 
treatment, we choose 

V a = V a0 + V a i , (f a = </9 a0 + (fi a l , 

(4.30) 

k B\ a = Bi a0 + Bi a i , B 2a = B 2a l + B 2a 2 , 

where the last subscript denotes both the order of approximation and the exponent 
of the power of t as t — > (remember that B 2 has an extra factor t as compared to 
Bi, cf. (2.11) (2.12)). The lowest approximation is defined by 

' id t V a0 + (l/2)AUaO = 

d t <Pajo = *"V(uao) + B 1L (v a0 ) 
< _ ; (4-31) 

BlaO = Bi{v a o) 

B 2a l = B 2 (V a0 , B la0 + Vifao) 

with initial conditions v a o(0) = v + , V'oo(l) = 0, so that in particular v a Q = U(t)v + 
and <f a0 = 0(1 — in t). Similarly the second approximation is defined by 

idtVai = remaining zero order terms of (2.16) not containing B 
d t (fai = 2t- l g{v a Q,v al ) + 2B 1L (v a0 , v al ) 

(4.32) 

Bioi = 2S 1 (v a0 ,Wal) 

B 2a 2 = order one terms of B 2 (v a , B la + B 2a Vip a ) 

with initial conditions v a i(0) = 0, (p a i(0) = 0. It turns out that the previous choice 
yields v a \ = 0(t(l — in t) 2 ), (p al = 0(t(l — in t) 2 ), and suffices to control the B 
independent part of the remainders. 

We now turn to the terms containing B , which are not covered by the previous 
choice. The difficulty comes from the fact that applying a derivative to B , whether 
in space or in time, generates a factor t~ l , as can be seen on the change of variable 
(2.6). Thus the most dangerous terms come from Vd t R\ which contains typically 
(Vd t B ) ■ Vi>+. Taking again h(t) = t 3 / 2+A for orientation, one needs an estimate 

II VdtR, \\ 2 <0(t- 1+x ) 
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whereas the best estimates on B for A a solution of the free wave equation yield 
only 

|| (Vd t B ) ■ Vv + || 2 < || Vd t B || 2 || Vv + iu < o (r 3 / 2 ) 



which is again off by slightly more than half a power of t. In order to deal with that 
problem, we impose a support condition on v + , namely 

Supp v + C {x : | \x\ — 1| > T]} (4.33) 

for some rj > 0. The effect of that condition is best understood when combined with 
B coming from a solution A of the free wave equation with compactly supported 
data. In fact, if 

Supp (A+,A+) C {x : |x| < R} , 
then by the Huyghens principle, A Q defined by (2.2) satisfies 

Supp A C {(x, t) : | \x\ - t\ < R} 

so that by (2.6) 

Supp B C {(x,t) : | \x\ — 1\ <t R} 
and therefore for all integers j, k, £, 

(V k diB ) VV = for t < T]/R ■ 

More generally, the support condition (4.33) together with some decay of (A + ,A + ) 
for large \x\ suffices to ensure the required estimates. We shall make the following 
assumption. 

(A5) Space decay of (A + ,A + ). Let \R De t ne characteristic function of the set 
{x : \x\ > R). We assume that (A + ,A + ) satisfies 

' || XR V k (x- V)SA+ \\ 2 V || XR V k (x- wy x - A + || 2 < c r- 1 

< 

k || xr{x- V) J 'i+;L 2 nL 6 / 5 || V || X r{x- V) J x-i + ;L 2 nL 6 / 5 || < C R' 1 

(4.34) 

for < j, k < 1 and for all R > Ro for some Rq > 0. 

The support condition (4.33) appeared in the early work [TT] on the MS system 
and was subsequently eliminated j2] [10] in the framework of the method of pj] . It 
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is an open question whether that condition can also be eliminated in the framework 
of the more complicated method described in the present paper. 
We can now state the result of the completion of Step 2. 

Proposition 4.2. Let v + G H 5 , xv + G H 4 and let v + satisfy the support condition 
(4-33). Let (A + ,A + ) satisfy the regularity and decay assumptions (A4) (A5). Let 
(v a ,ip a ,B a ) be defined by (4.30)-(4.32) with v a0 (0) = v+, ip a0 {l) = 0, v al (0) = 0, 
</ 3 oi(0) = and B 0a = B . Then the assumptions (Al)-(AS) are satisfied with 

h{t) =t 2 (l- Int) 4 . (4.35) 

Putting together Propositions 4.1 and 4.2 yields the main result on the Cauchy 
problem at t — for the auxiliary system (2.16)-(2.18). 

Proposition 4.3. Let v + G H 5 , xv + G H 4 and let v + satisfy the support condition 
(4-33). Let (A + ,A + ) satisfy the regularity and decay assumptions (A4) (A5). Let 
(v a ,¥a,B a ) be defined by (4.30)-(4.32) with v a0 (0) = v + , <p a0 (l) = 0, v al (0) = 0, 
<^oi(0) = and B 0a = B . Let h(t) = t 2 (l — in t) 4 . Then there exists r, < r < I, 
and there exists a unique solution (v , <p, B 2 ) of the auxiliary system (2.16)-(2.18) 
such that (v — v a , B 2 — B 2a ) G -X"((0, r]) and such that (<p — (p a )(t) ^0 as t — > 0. 

We now return to the Cauchy problem at infinity for the MS system in the form 
(1.3) (2.1). We start from the asymptotic data (w + , A + , A + ) for (u, A) and we define 
v + = Fu^. We define A by (2.2) and B by (2.6) . We define (v a , ip a , B a ) by (4.30)- 
(4.32) with the appropriate initial conditions. We then define the asymptotic form 
(u a ,A a ) of (u, A) in analogy with (2.4)-(2.6), (2.14) by 

u ac = v a exp (-itp a ) , (4.36) 

u a (t) = Fu ac {l/t) , (4.37) 

A a (t) = -r l D (t) B a (l/t) . (4.38) 

The final result can then be stated as follows. We recall that u is defined by 
(2.5), that V* is defined by (3.2), that S is defined by (3.3) and that j, k, I denote 
non negative integers. 
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Proposition 4.4. Let u + be such that v + = Fu + G -ff 5 rcn + G id 4 and £na£ n + 
satisfy the support condition (4-33). Let (A + ,A + ) satisfy the regularity and decay 
assumptions (A4) (A5). Let (u a ,A a ) be defined by (4-36)-(4-38) with (v a ,tp a ,B a ) 
defined by (4.30)-(4.32) with v a0 (0) = v + , </? a0 (l) = 0, n a i(0) = 0, y? a i(0) = and 
Bo a = Bq. Then there exists T > 1 and there exists a unique solution (u,A) of the 
MS system (1.3) (2.1) in I = [T, oo) with the following properties : 

u G C(I, V*) , d t u G C(I, H 1 n FH 1 ) , 

A G C(I, H 1 n H 2 ) , x-Ae C(I, H 2 ) , SA G C(I, H 1 ) , 
and (u, A) satisfy the estimates 

II w(f);K || < C(l + £n t) 3 (4.39) 

|| V(A-A o )(0 || 2 < Ct' 1 ' 2 (4.40) 

/or all t G /. Furthermore (u, A) behaves asymptotically as (u a , A a ) in the sense 
that the following estimates hold for all t & I : 

|| x k V\u - u a )(t) \\ 2 <C r 2+k/2 (l + £n tf (4.41) 

fori = 0, 1 and < fc + i < 3, 

|| ^ t (n-5 a )(t) || 2 < C r 3+k / 2 (l+£nt) 4 (4.42) 

/or k = 0,1, 

|| V fe+1 ^'(A - A a )(t) || 2 V r 1 1| V k+l &x ■ (A - Aa)(t) \\ 2 

< c r 5/2-fc/2( 1+ ^ n f )4 

/or < j, /c, j + A; < 1. 



The existence part of the proof of Proposition 4.4 is obtained by first constructing 
(v, (p, B 2 ) by the use of Proposition 4.3, then constructing (n, A) from (v, ip, B 2 ) by 
using the change of variables (2.4)-(2.6) (2.10) (2.11) (2.14), and finally translating 
the properties of (v, ip, B 2 ) that follow from Proposition 4.3 in terms of (n, A). The 
uniqueness part of the proof follows readily from Proposition 3.2. In fact, if (tij, Ai), 
% — 1,2, are two solutions of the system (1.3) (2.1) obtained from Proposition 4.4 
with the same (u a ,A a ), then the conditions (3.9) (3.10) follow from (4.39) (4.40) 
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with a = 3, while the condition (4.41) together with the commutation relation 
xU(-t) = U(-t)(x + itV) implies 

||< x/t > (uj - u 2 )(t) || 2 < C r 2 (l + in tf 

which implies (3.11). 
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